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Abstract. A large set of TREESPH simulations is used to test the global 
morphology of galaxy clusters and its evolution against X-ray data. A powerful 
method to investigate substructures in galaxy clusters are the power ratios in- 
troduced by Buote & Tsai. We consider three flat cosmological models: CDM, 
ACDM (f^A = 0.7) and CHDM (flh = 0.2, 1 massive ly), all normalized so to 
fit the observed number of clusters. For each model we built 40 clusters, using 
a TREESPH code, and performed a statistical comparison with a data sample 
including nearby clusters observed with ROSAT PSPC instrument. The com- 
parison disfavors the ACDM model, as clusters appear too relaxed, while CDM 
and CHDM clusters, in which a higher degree of complexity occurs, seem to be 
closer to observations. A better fit of data can be expected for some different DM 
mix. If DM distributions are used instead of baryons, we find substructures more 
pronounced than in gas and models have a different score. Using hydrodynamical 
simulations is therefore essential to our aims. 



1 Introduction 

Clusters of galaxies are likely to be dynamically young systems and a promising 
way to constrain cosmological models arises from the study their substructures 
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Among the methods suggested to quantify the degree of inhomogeneity in 
clusters R; Q ), one of the most promising is based on the so-called power 

ratios n^^' ([l^], hereafter VGB99; |^). It amounts to a multipole expansion 
accounting for the angle dependence of cluster X-ray surface brightness, limited 
to the first few multipole terms and at fixed scales (typically of the order of the 
Mpc). This is an effective and synthetic way to discriminate cluster features and 
H^"'-' are found do depend on the cosmological model enabling to discriminate 
among different cosmologies. 



2 Power ratios: definition and evolution 



The number of photons coUected by ROSAT PSPC does not depend on the tem- 
perature r, if it is ^ 1 keV. We can then assume that the X-ray surface brightness 
Ex ~ A J pl{r)dz. Here A ~ const results from an integration along the line of 
sight. 

The procedure to work out If'™' is as follows (see, e.g., 1^, hereafter G98; 
VGB99; ): (i) Pt{r) is projected along a line of sight on a (random) plane to yield 
the X-ray surface brightness E(i?, Lp); the centroid is used as origin, (ii) By solving 
the Poisson equation V^"l? = E(_R, ip) we obtain the pseudo-potential $(-R, (p). (iii) 
The coefficients of the expansion of $ in plane harmonics will be used to build the 
power ratios n'-"'\Rap) = logio(-Pm/fb). Here, Pm{R) = {al,+(5l)/2m^ ,Po = 
[ao ln(i?/kpc)]^, while 

ctm = I ds s^^^ I dip[T,{sR,ip)R^]cos{mp), (1) 
Jo Jo 

and (3m. has an identical definition, with sin instead of cos. Owing to the definition 
of the centroid, H'^' vanishes. We shall restrict our analysis to n''"'(m = 2, 3, 4), 
to account for substructures on scales not much below R itself. We consider three 
different aperture radii Rap = 0.4,0.8, i.2/i^^Mpc. 

Because of its evolution, a cluster moves along a curve of the 3-dimensional 
space spanned by such fl^'^^'s; this curve is called evolutionary track. Quite in 
general, a cluster starts from a configuration away from the origin, corresponding 
to a large amount of internal structure and evolves towards isotropization and 
homogeneization. 

Actual data, of course, do not follow the motion of a given cluster along the 
evolutionary track. Different clusters, however, lie at different redshifts and can 
be used to describe a succession of evolutionary stages. 

3 The simulated and the observed cluster sample 

We consider three spatially flat cosmological models: CDM, ACDM with a cos- 
mological constant accounting for 70% of the critical density, and CHDM with 
1 massive neutrino with mass = 4.65 eV, yielding a HDM density parameter 

= 0.20. We set h = 0.5 for CDM and CHDM and h = 0.7 for ACDM; for all 
models the primeval spectral index n — 1 and the baryon density parameter is 
selected to give Qtfi^ = 0.015. All models were normalized in order to reproduce 
the present observed cluster abundance ([^, [^). 

In order to achieve a safe statistical basis for our analysis, for each cosmological 
model, we select the 40 most massive clusters from an N-body P3M simulation. 
For each of them, we perform a hydrodynamical TREESPH simulation (VGB99, 
G98; see also § and |ll)). 
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Figure 1: Histograms of iP in the Rap = 0.8h ^Mpc for simulated clusters. 

Clusters are distributed in redshift so to reproduce the same redshfit distribu- 
tion of the observed cluster sample. 

Our observed data set is the same used by including nearby {z ^ 0.2) 
clusters observed with ROSAT PSPC instrument (see VGB99, G98 for details). 
The resulting sample is partially incomplete, but, clusters were not selected for 
reasons related to their morphology and the missing clusters are expected to have 
a distribution of power ratios similar to the observed one. 

4 Results and Conclusions 

For simulated clusters, power ratios II'™'-' have been computed from the gas dis- 
tribution. 

A visual inspection of how the n'™' are distributed can be obtained from 
Fig.l, whose histograms show the fraction of clusters with a given n^™' . Here we 
show the distribution for 11'' and Rap = 0.8/i~^Mpc for each cosmological model 
and for the ROSAT data sample; distributions for the other n^*"^ and the other 
apertures show a similar behavior 

Quite in general we can conclude that while CDM and CHDM are marginally 
consistent with data, ACDM is far below them. In order to quantify these dif- 
ferences, we used the Student t-test, the F-test and the Kolgomorov-Smirnov 
(KS) test. For example, according to the t-test, the probability p-t that the 
simulated and observed power ratios distributions are originated from the same 
process is roughly in the ranges 0.11-0.60 for CDM, 0.03-0.37 for CHDMand 
0.15-10~*-0.87-10~^ for ACDM. The other statistical tests provide similar proba- 
bilities. Such figures seem to exclude that ACDM can be considered a reasonable 
approximation to data. The best score belongs to CDM, but also CHDM is not 
fully excluded and different mixtures could certainly have better performance. 
An inspection of the model clusters actually shows that the ACDM model does 
produce less substructures than the other models do. A possible interpretation of 
such output is that the actual amount of substructures is governed by Oq rather 




than by the shape of power spectra. 

According to the same tests, if cosinological models are compared with data 
on the basis of DM n'"'^ , values are shifted, indicating an increase in the amount 
of substructures for DM with respect to the gas. This is to be ascribed to the 
smoothing effects of the interactions among gas particles, which erase anisotropics 
and structures, while DM Il'^^scarcely feel dissipative processes. Hence, using 
DM n^'^Heads to biased scores: CDM and CHDM models keep too many sub- 
structures and axe no longer consistent with data; on the contrary, the increase of 
substructures pushes ACDM to agree with ROSAT sample outputs. 

We also considered the cluster distribution in the 3 dimensional parameter 
space with a^xes given by II'™' (m = 2, 3, 4), as well as projections of such distri- 
butions on planes. Comparing such distributions for data and models, we find a 
significantly stronger correlation of n'™' in models than in data. Distributions for 
simulated clusters show a linear trend while distributions of observed clusters tend 
to be more scattered than simulated points. The degree of correlation depends on 
the model, but seems however in disagreement with data. Model clusters tend to 
indicate a significantly faster evolution than data. The cosmological model which 
seems closest to data is CHDM and it is possible that different CHDM mixtures 
can lead to further improvements. Also A models with flrn > 0.5 might deserve 
to be explored. Virializcd clusters had their turn-around at a time < to/3 {to'. 
present age of the Universe). In turn, Qa becomes dominant at 2 = (fiA/f^m)'^'^^- 
If such redshift occurs at a time ~ to/3, we expect results from A models to be 
closer to observations. 
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